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Abstract 


Data  archives  are  such  that,  for  a  selected  cruising  level,  a  reasonably  large 
sample  of  daily  average  temperatures  along  a  specific  route  can  be  obtained  and 
summarized  into  a  frequency  distribution.  But  if,  as  in  Air  Force  operations, 
many  different  and  constantly  changing  routes  are  involved,  then  the  climatologic 
tools  and  methods  must  be  such  as  to  provide  estimates  of  the  route  average  temper 
ature  without  recourse  to  the  archives  and  the  processing  of  raw  data.  A  method 
has  been  devised  that  utilizes  a  model  of  the  horizontal  spatial  correlations  of 
temperature.  It  has  been  demonstrated  on  a  continent -wide  area,  and  shows 
promise  of  hemisphere -wide  application. 
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A  Statistical  Model  of  the  Temperature  Field  at 

Supersonic  Aircraft  Altitudes 


I.  INTROIJUCTION 

Supersonic  aircraft  (SSA)  will  cruise  at  60000-  to  80000-ft  altitudes  where 
the  ambient  air  temperature  will  affect  the  rate  of  fuel  consumption.  It  has  been 
estimated  (Nelms,  1964),  with  some  qualification  (Stickle,  1965),  that  the  allowable 
payload  on  a  3000-mile  trip  will  be  reduced  by  perhaps  nine  passengers  and  their 
baggage  when  the  average  route  temperature  at  cruising  level  is  20°F  warmer 
than  the  U.  S.  Standard.  At  the  same  time  it  appears  that  a  flying  mission  will 
not  gain  by  a  diversion  from  a  great -circle  route  to  avoid  extremes  of  route 
temperature  (Serebreny,  1964). 

The  subject  of  the  average,  or  integrated,  route  temperature  has  not  been 
neglected  in  the  past  half-dozen  years  (Crutcher,  1963;  Charles,  1964;  Court  and 
Abrahms,  1964).  The  U.  S.  Weather  Bureau,  under  FAA  sponsorship,  issued  three 
special  reports  on  the  routes  New  York  to  San  Francisco  (1963),  San  Francisco  to 
Thule  to  Stockholm  (1964)  and  New  York  to  Paris  (1964),  in  which  there  are  tables 
of  the  distribution  of  average  route  temperatures  at  several  key  levels  for  all 
seasons.  The  method  that  the  USWB  used  to  obtain  the  route  averages  was  unequi¬ 
vocally  sound.  On  each  of  a  large  sample  of  days,  at  each  level  and  season,  the 
average  temperature  along  the  selected  route  was  computed.  In  Figure  1,  curve 
(c)  shows  the  frequency  distribution  of  these  averages  for  the  San  Franciso  to  Thule 
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TEMPERATURE  T  (°C)-* 

Figure  1.  The  Cumulative  Frequency  Distribution  P(T)  of  Temperature  T, 

50  mb.  Winter,  Over  (a)  San  Francisco,  and  (b)  Thule;  (c)  San  Francisco  — 
Thule  Route  Average 

> 

route  in  winter.  Boeing  Aircraft  Co.  , 'Wells,  1962  a,  b)  has  issued  two  reports 
on  the  average  route  temperature  meins  and  standard  deviations  for  many  conceiv¬ 
able  routes  throughout  the  world.  All  these  tables,  however,  have  to  date  been 
limited  to  the  levels  of  subsonic  jet  aircraft.  Neal  Barr,  of  Boeing,  has  in  a 
private  communication  mentioned  that  a  revision  that  will  include  the  50-mb  level 
is  about  to  be  published. 

From  the  standpoint  of  Air  Force  operations,  consideration  must  not  be 
restricted  to  a  set  of  preselected  routes.  It  must  be  possible  to  find  the  probability 
of  temperature  —  warm  or  cold  —  at  relatively  short  notice  for  any  route  and 
without  recourse  to  the  data  archives  and  a  computer  effort.  It  has  therefore  been 
deemed  necessary  to  devise  a  model  whereby  temperatures  can  be  related  spatially. 
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On  the  assumption  that  the  SSA  will  fly  at  a  constant  pressure  height,  the  model 
described  below  has  been  applied  to  the  North  American  continent,  at  levels  50 -mb 
and  30-mb  (67500  and  78000  ft)  in  all  four  seasons  of  the  year. 


2.  THE  BASIC  RELATION 

At  any  given  altitude  and  season,  the  temperature  frequency  distribution  over 
a  given  station  generally  departs  from  a  normal  or  gaussian  distribution  (Barr,  1964). 
More  often  than  not  this  departure  is  small,  but  sometimes  it  becomes  large, 
especially  over  the  polar  regions  at  altitudes  of  60000  to  80000  ft  where  the 
temperature  distribution  in  winter  is  bimodal,  owing  to  a  well -documented 
phenomenon  termed  the  explosive  warming  (Reed  et  al. ,  1963;  Craig  and  Hering, 

1959).  When  the  temperature  is  averaged  over  a  route  of  2000  naut.  mLor  longer, 
however,  a  centralizing  effect  produces  a  more  nearly  normal  distribution.  For 
example.  Figure  1  shows  the  cumulative  frequency  distribution  of  the  temperature 
in  winter  at  50  mb:  (a)  over  San  Francisco,  (b)  over  Thule,  Greenland.  Curve  (b) 
shows  distinct  bimodality,  but  curve  (c),the  distribution  of  the  average  over  the 
entire  route  from  San  Francisco  to  Thule,  is  much  less  bimodal.  For  the  model 
described  below,  all  route -integrated  temperatures  are  assumed  to  be  normally 
distributed  and  therefore  adequately  defined  by  the  mean  and  the  standard  deviation 
of  the  distribution. 

Consider  a  route  defined  by  (ml)  checkpoints  Xq,  ...  ,Xn,  including  the  point  of 
departure  XQ  and  the  destination  Xn.  Let  sq,  ...  ,  sn>  be  the  weights  to  be  attached 
to  each  observation  for  a  suitable  determination  of  route  average  temperature,  and 
T0. ...» Tn,  the  temperature  of  each  of  the  (n+1)  checkpoints.  Then  the  route  average 
T  is  given  by 

n  n 

T  5  £  eiTi/]C  v  <l) 

i=o  i=o 

to  give  the  mean 

t  ■  e  si  vix  <2> 

i=o  i=o 

The  standard  deviation  a  of  the  route  average  is  given  in  terms  of  the  standard 
deviation  of  each  temperature  T^  and  the  correlation  between  the  ith  and 
jth  terms  in  the  average.  Thus, 


b 
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*2  ‘ 

i,j=0  '  '  i=o 

When  a  route  is  divided  into  relatively  short  segments  by  the  checkpoints  Xq,  ...  ,Xn, 
the  weight  sQ  should  be  taken  as  the  length  from  XQ  to  midway  between  XQ  and 
Xj  ,  the  weight  should  be  taken  as  the  length  from  midway  between  XQ  and  Xj 
to  midway  between  X  ^  and  Xg.  and  so  on.  Table  1  lists  the  winter  50-mb  mean  Tj 
and  standard  deviation  oj  at  each  of  eight  checkpoints  on  the  route  from  San  Fran¬ 
cisco  to  Thule.  The  total  distance  of  2810  naut.  mi.  was  divided  into  eight  lengths, 
shown  as  the  weighting  factors  s^.  The  data  sample  that  yielded  the  means  and 
standard  deviations  (1958  to  1963)  also  yielded  the  correlation  coefficients  in  the 
table.  Thus,  Table  1  provides  the  information  necessary  and  sufficient  for  esti¬ 
mating  mean  T  and  standard  deviation  a  by  means  of  Kqs.  (2)  and  (3). 


Table  1.  The  San  Francisco  — Thule  Route  Temperatures,  Winter,  50  mb,  at 
Eight  Checkpoints.  (Tj  s  temperature  mean,  in  °C;  -  standard  deviation,  in 

°Cj  sj  =  length  of  the  route,  in  nautical  miles,  represented  by  the  temperature 
at  the  checkpoints.  The  body  of  the  table  contains  correlation  coefficients  between 
temperatures  of  pairs  of  checkpoints.  ) 


Station  No. 
Location 

0 

SF 

1 

Medford 

2 

Seattle 

3 

Edmonton 

4 

Ft.  Smith 

5 

Baker  Lake 

6 

Resolute 

7  1 

Thule 

Mean  T^ 

-60 

-59 

-58 

-57 

-57 

-57 

-62 

-62 

SD  a. 

3.  3 

3.  8 

4.  1 

6.  0 

8.  4 

8.  8 

10.  6 

9.  7 

Length  s. 

150 

320 

390 

390 

400 

480 

470 

210 

Station  No. 

0 

1. 000 

.  85 

.  73 

.  43 

.  24 

.  05 

-.  09 

-.  27 

1 

1.  000 

.  89 

.  63 

.  53 

.  19 

.  00 

-.  28 

2 

1.  000 

.  80 

.  68 

.  22 

.  12 

-.  26 

3 

1.  000 

.  90 

.  66 

.  40 

.  06 

4 

1. 000 

.  89 

.  69 

.  38 

|  5 

1.  000 

.  85 

.  62 

1  6 

1.  000 

.  89 

7 

1. 000 

This  technique  yielded  the  means  of  the  average  route  temperatures  and 
estimated  standard  deviations  in  Column  B,  Table  2.  For  comparison,  the  FAA 
standard  deviation  values  based  on  1957  to  1962  data,  computed  directly  from  day- 
to-day  route  averages,  are  shown  in  Column  C. 


Both  Column  B  and  C  estimates  are  made  without  modeling  and  are  equally 
valid.  The  difference  is  therefore  an  indication  of  the  error  in  estimating  the 
route -temperature  standard  deviation.  The  root -mean -square  difference  of  the 
four  pairs  of  values  is  0.  4°C. 

The  same  technique  could  be  applied  to  obtain  the  mean  and  standard  deviation 
of  the  route  average  for  any  hemispheric  route.  Success  would  depend  on  using  a 
sufficient  number  of  route  checkpoints  and  knowing  not  only  the  individual  means 
and  standard  deviations  of  each  checkpoint  but  the  correlation  coefficients  between 
pairs.  The  correlation  coefficients  present  a  difficulty.  To  adequately  cover  the 
Northern  Hemisphere,  or  even  the  North  American  quadrant,  the  tables  to  give 
such  intercorrelations  would  have  to  be  voluminous,  to  say  nothing  of  the  massive 
effort  that  would  be  required  to  produce  them.  It  has  therefore  been  found  necessary 
to  develop  a  model  that  yields  usable  estimates  of  these  correlation  coefficients. 


Table  2,  Comparison  of  Means  and  Standard  Deviations 
of  Rqute  Average  Temperatures  at  50  mb 


Season 

Route 

Mean 

<°C) 

Standard  Deviation  (°C) 

A* 

B+ 

CT 

Summer 

SF  -  NY 

-53.  5 

.  98 

.  96 

1.  00 

SF  -  Thule 

-48 

1.  22 

1.  26 

1.  17 

W  inter 

SF  -  NY 

-59 

2.  70 

2.  62 

2.  36 

SF  -  Thule 

-59 

5.  63 

5.  68 

5.  48 

*by  model 

"hjy  formula  (3) 

+from  USWB  reports 


3.  THE  MODI'] I . 

Figure  2  is  an  example  of  the  kind  of  field  that  a  model  must  be  able  to  approx - 

mate.  The  lines  join  points  having  equal  correlation  coefficients,  based  on  the 

50-mb  winter  temperature  at  the  master  station,  San  Francisco,  and  the  continental 

slave  stations.  As  might  be  expected,  the  correlation  close  to  the  master  station 

is  high.  It  decreases  to  zero  at  about  1800  naut.  mi.  from  the  master  station,  and 

becomes  slightly  negative  at  greater  distances. 

After  some  trial  and  error,  a  model  with  a  single  parameter  (o  y)  was  chosen 

to  estimate  the  correlation  p..  between  points  X.  and  X  .  .  If  the  distance  between 

1J  ^  J 

these  points  is  y^  as  measured  in  earth  radians  on  a  polar  stereographic  chart, 
then  the  estimate  of  is 
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Figure  2.  Isopleths  of  Correlation  Coefficient  Between  the  50 -mb 
Temperature,  Winter,  Over  San  Francisco  and  Other  North  Ameri¬ 
can  Stations  f  projection:  polar  stereographic;  scale  measure:  earth 
radians  (tenths)] 
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exp 


coa(aij  V 


as  drawn  in  Figure  3. 


(4) 


Figure  3.  Graph  of  the  Model  of  Correlation  Coefficient  (p. .) 
Decay  With  Distance  [yy  *  Distance  in  Radians  Between  ith  J 
and  jth  Stations j  u  =  Model  Parameter] 


The  model  in  Figure  3  [Eq.  (4)],  has  merit  but  it  also  has  limitations.  The 
correlation  coefficient  will  decrease  to  zero  and  become  somewhat  negative;  it  cannot 
be  more  negative  than  -0.  06,  but  this  is  not  an  important  consideration  when  the  cHstance 
between  stations  is  great.  Of  greater  significance,  the  more  independent  the  temp  rature 
over  a  given  station,  the  larger  the  parameter  a.  j  .  Since  the  distance  y^  is  mo*  ie  equal 
to  the  distance  as  measured  on  a  polar  stereographic  projection  instead  of  equal  to 
the  actual  terrestrial  distance  between  stations,  the  effect  of  the  varying  scale  is 
absorbed  in  the  parameter  a.j  . 

To  construct  the  model,  let  the  origin  on  a  polar  stereographic  chart  be  at  the 

N  pole,  the  x  axis  positive  along  the  0°  meridian,  and  the  y  axis  positive  along 

the  90°E  meridian.  Let  X  denote  longitude  and  <f>  denote  latitude.  Then  the 

distance  y^  between  the  two  points  (x.,  y.)  and  (x  ,  y  ),  which  respectively 
J  J  J 

represent  the  two  points  and  (A.,0.),  is  given  by 

J  J 


(5) 


ij 


■V 


(xj 


(yj  -y^ 


where 


<p. 

=  2  cos  A.  tan(|-  -  , 

_  <p: 

yi  =  2  sin  A.  tan  (|-  - 
Likewise, 


Xj  =  2  cos  A  j  tan(|-  -  ^) , 
y-  =  2  sin  Aj  tan(|-  -  tp)  . 

Thus,  the  radial  distance  from  pole  to  equator,  which  is  t/2  radians  on  the 
earth's  surface,  is  2.  0  radians  on  the  polar  stereographic  chart.  Figure  2  contains 
a  constant  scale  for  the  chart,  marked  off  in  tenths  of  radians  (see  Table  3). 


Table  3.  Construction  of  Scale  Yielding  Segments 
«  of  Equal  Length  on  the  Polar  Stereographic  Chart 
[All  measurements  are  made  from  the  N  pole  along 
the  0°  meridian  (  x  axis).  ] 


y  =  o;  x  = 

2  tnn(£  -  £) 

X 

</.(°N) 

0.  1 

84.  2 

0.  2 

78.  G 

0.  3 

G7.  4 

j  0.  4 

G7.  4 

0.  5 

G2.  0 

0.  G 

5G.  6 

0.  7 

51.  4 

0.  8 

4G.  4 

0.  9 

41.  G 

1.  o 

30.  8 
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At  this  stage  the  model  does  not  appear  especially  useful.  The  parameter  a- 
is  a  function  of  the  two  stations  X^  and  Xj  .  What  is  needed  is  a  single  parameter 
that  can  be  assigned  to  each  station  individually  so  that  it  can  be  mapped  in  much 
the  same  way  as  the  mean  or  standard  deviation  is  mapped. 

Arbitrarily  assign  the  number  to  the  station  and  make  the  parameter 
the  average  of  two  values: 


(ai  +  cr  ■) 


(6) 


But  this  arbitrary  settlement  of  the  question  leads  to  another  question:  What  is  the 
best  value  of  a,  to  assign  to  station  Xi  ?  Clearly,  it  should  be  a  value  that  mini¬ 
mizes  the  error  in  estimation  of  p  by  Eqs.  (4)  and  (6). 

To  attack  the  problem  expeditiously  it  was  decided  to  find  the  value  of  at 
each  master  station  X^  so  that  when  was  used  in  place  of  in  Eq.  (4),  the 

result  would  be  a  set  of  estimates  of  p  between  the  ith  station  and  all  other 

stations  within  a  reasonable  distance  of  the  ith  station.  The  value  of  a.  was 
chosen  such  that 


E. 

x 


65 

z 

j=l 


w 


.. (p..  -$.y 


(7) 


was  a  minimum,  where  E  is  the  weighted  sum  of  squares  of  error.  The  quantity 
Wy  is  a  weighting  factor.  It  can  be  seen  from  Eq.  (3)  and  its  application  in  Table  1 
that  it  is  more  important  to  closely  estimate  the  correlations  of  neighboring  stations 
than  to  estimate  the  correlations  of  stations  farther  apart.  After  a  few  trials,  it 
was  decided  to  set  the  weighting  factor 


=  0.6  -yy.  yJJ<0.6. 

=  0  ,  y .j  >  0.  6 . 


(8) 


4.  DATA  AM)  RESULTS 

To  apply  the  model,  the  semidaily  temperatures  for  the  five  years  1959  to  1963 
were  assembled  for  the  50-  and  30-mb  heights  at  which  SST  are  expected  to  cruise. 
The  data,  from  66  stations  over  North  and  Central  America  (Table  4),  were  grouped 


Table  4.  The  Sixty -six  Stations  Whose  Observations  of  the  50-mb 
and  30-mb  Temperatures  for  the  Period  of  Record,  December  1958 
to  November  1963,  Were  Correlated  forThis  Report 


Sta. 

No. 

Location 

Latitude 

Longitude 

1 

San  Diego,  California 

32 

49 

117 

08 

2 

Balboa,  Panama,  Canal  Zone 

8 

59 

079 

34 

3 

San  Juan,  Puerto  Rico 

18 

26 

066 

00 

4 

Grand  Cayman,  BWI/Roberts  Fid 

19 

18 

081 

22 

5 

San  Andrei,  Colombia 

12 

35 

081 

42 

6 

Miami,  Florida 

25 

48 

080 

16 

7 

Burrwood,  Louisiana 

28 

58 

089 

22 

8 

Merida,  Mexico 

20 

55 

089 

40 

9 

Brownsville,  Texas 

25 

54 

097 

26 

10 

Greensboro,  North  Carolina 

36 

05 

071 

57 

11 

Charleston,  South  Carolina 

32 

54 

080 

02 

12 

Nashville,  Tennessee 

36 

07 

086 

41 

13 

Fort  Worth,  Texas  /Carswell  AFB 

32 

50 

097 

03 

14 

Jackson,  Mississippi 

32 

20 

090 

13 

15 

Dodge  City,  Kansas 

37 

46 

099 

58 

16 

Sable  Island,  Nova  Scotia 

43 

56 

060 

02 

17 

Nantucket,  Massachusetts 

41 

15 

070 

04 

18 

Portland,  Maine 

43 

39 

070 

19 

19 

Peoria,  Illinois 

40 

40 

089 

41 

20 

Sault  Ste  Marie,  Michigan 

46 

28 

084 

22 

21 

Goose,  Newfoundland 

53 

18 

060 

27 

22 

Fort  Chimo,  Quebec 

58 

06 

068 

26 

23 

Nitchequon,  Quebec 

53 

12 

070 

54 

24 

Fort  Harrison,  Quebec 

58 

27 

078 

08 

25 

Moosonee,  Ontario 

51 

16 

080 

39 

26 

Trout  Lake,  Ontario 

53 

50 

089 

52 

27 

Churchill,  Manitoba 

58 

45 

094 

04 

28 

Frobisher  Bay,  Northwest  Territory 

63 

45 

068 

33 

29 

Hall  Lake,  Northwest  Territory 

68 

47 

087 

15 

30 

Baker  Lake,  Northwest  Territory 

64 

18 

096 

00 

31 

Thule,  Greenland 

76 

33 

068 

49 

32 

Resolute  Bay,  Northwest  Territory 

74 

43 

094 

59 

33 

Alert,  Northwest  Territory 

82 

30 

062 

20 

34 

Eureka  Sound,  Northwest  Territory 

80 

00 

085 

50 

35 

Del  Rio,  Texas 

29 

23 

100 

46 

36 

Midland,  Texas/Sloan  Fid 

32 

56 

102 

12 

37 

El  Paso,  Texas 

31 

48 

106 

24 

38 

Grand  Junction,  Colorado 

39 

07 

108 

32 

39 

Ely,  Nevada 

39 

17 

114 

51 

40 

Tucson,  Arizona 

32 

07 

110 

56 

41 

Oakland,  California 

37 

44 

122 

12 

42 

Bismarck,  North  Dakota 

46 

46 

100 

45 

43 

North  Platte,  Nebraska 

41 

08 

100 

41 

44 

Great  Falls,  Montana 

47 

29 

111 

21 

45 

Medford,  Oregon 

42 

22 

122 

52 

46 

Olympia,  Washington 

47 

27 

122 

18 

47 

The  Pas,  Manitoba 

53 

58 

101 

06 

48 

Edmonton,  Alberta 

53 

34 

113 

31 

49 

Prince  George,  British  Columbia  I 

53 

50 

122 

41 

50 

Fort  '.'elson,  British  Columbia 

58 

50 

122 

35 

51 

Annette,  Alaska 

55 

02 

131 

34 

52 

Yakutat,  Alaska 

59 

31 

139 

40 

53 

Cold  Bay,  Alaska 

55 

12 

162 

43 

54 

St  Paul  Island,  Alaska 

57 

09 

170 

13 

55 

Fort  Smith,  Northwest  Territory 

60 

01 

111 

58 

56 

Coppermine,  Northwest  Territory 

67 

49 

115 

05 

57 

Anchorage,  Alaska 

61 

10 

149 

59 

58 

Fairbanks,  Alaska 

64 

49 

147 

52 

59 

Bethel,  Alaska 

60 

47 

161 

48 

60 

Nome,  Alaska 

64 

30 

165 

26 

61 

Isachsen,  Northwest  Territory 

78 

47 

103 

32 

62 

Mould  Bay,  Northwest  Territory 

76 

14 

119 

20 

63 

Barter  Island,  Alaska 

70 

08 

143 

38 

64 

Barrow,  Alaska 

71 

18 

156 

07 

65 

Shemya,  Alaska 

52 

43 

174 

06 

66 

Pittsburgh,  Pennsylvania 

40 

30 

080 

13 

by  seasons  (winter:  December,  January,  February;  summer:  June,  July,  August). 
For  each  station  in  each  season  the  sample  mean  and  standard  deviation  of  the 
temperature  were  found  and  the  correlation  coefficient  between  pairs  of  all  stations 
calculated.  The  means  and  standard  deviations  were  then  plotted  and  isopleths 
drawn  (Figures  4  to  11,  a,  b).  Equation  (4)  was  used  in  making  the  estimate 
(a^)  and  the  error  of  estimate 


was  squared,  multiplied  by  and  added  with  other  squared  errors  to  give,  for 

the  ith  station: 


By  trying  a  set  of  values  for  a  .  , 
a .  -  0.  5,  1.5,...,  10.  5  , 

it  was  possible  to  find  the  value  of  that  minimized  E,,.  The  results  for  a  are 
shown  in  Figures  4  to  11  (c). 

In  determining  the  sample  mean  T.,  standard  deviation  u.,  and  the  correlation 
coefficient  p..»  between  stations,  the  diurnal  cycle  was  ignored.  To  find  out  how 
serious  an  oversight  this  could  be,  the  OOOOZ  and  1200Z  temperatures  were  averaged 
separately  and  their  difference  D(D  =  T120q  -  T0000)  plotted  for  the  winter  and 
summer  seasons  at  50  mb  (Figures  12  and  13).  The  effect  of  the  daytime  heating 
at  the  50 -mb  level  is  apparent.  Close  to  the  90°W  meridian  the  difference  is  zero; 
to  the  east  it  is  positive  and  to  the  west  it  is  negative.  It  is  not,  however,  clear 
whether  the  effect  is  real  or  due  to  heating  of  the  instrument.  This  subject  has 
been  dealt  with  in  some  detail  by  Chiu  ( 1959),  Badgley  (1959),  and  Harris  et  al. 

(1962). 

Between  the  standard  deviations  of  OOOOZ  and  1200Z  observations,  there  is  a 
root -mean-square  difference  of  0.  124°C  in  summer  and  0.  272°C  in  winter.  Errors 
of  observation  should  increase  the  observed  variances  by  a  mean-square  amount 

2  2  2 
e  a  (observed)  °  (true) 


The  value  of  e  has  been  found  to  be  ( Johannessen,  1959)  approximately  0.  7.  By 
comparison,  the  root -mean -square  difference  between  OOOOZ  and  1200Z  observa¬ 
tions  can  be  considered  negligible,  whether  it  is  real  or  due  to  errors  of  the 
instrument. 
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For  the  present  study,  since  the  standard  deviation  is  of  the  order  of  2°  to 
10°C,  all  errors  of  observation  are  considered  to  have  only  slight  effect  on  the 
values  of  the  standard  deviation  and  are  therefore  neglected.  Eventually,  refined 
observations  and  larger  samples  of  data  could  be  used  to  yield  improved  values  of 
means,  standard  deviations,  and  correlation  coefficients,  so  that  the  charts  of  this 
paper  might  be  improved.  There  are  similarities,  but  also  significant  differences, 
between  the  charts  in  this  report  and  the  charts  or  tables  of  Ratner  (1957),  Air 
Weather  Service  (1953),  Li  Peng  (1965),  and  Smith,  McMurray,  and  Crutcher  (1963). 
All  the  charts  differ  from  one  another  in  significant  values  but  the  50-mb  and  30-mb 
averages  given  by  Ratner  from  the  1946-1955  records  are  in  satisfactory  agreement 
with  Figures  4a  to  11a  in  this  report. 


5.  APPLICATION 

From  the  (a)  and  (b)  charts  of  Figures  4  to  11  we  read  the  mean  and  standard 
deviation  of  temperatures  at  checkpoints  along  a  selected  route.  From  the  (c) 
charts  of  Figures  4  to  11  we  read  the  value  of  a  at  each  checkpoint,  and  th  n 
obtain  the  correlation  coefficients  between  checkpoints.  This  provides  the  requisite 
information  for  obtaining  the  mean  and  standard  deviation  of  the  route  average.  The 
results  obtained  for  the  San  Francisco  —  Thule  route,  winter,  50  mb  (Table  5),  illus¬ 
trate  the  procedure. 

The  chart  route  is  divided  into  equal  segments  0.  1  radian  long  (Table  3).  For 
the  San  Francisco  —  Thule  route  this  calls  for  nine  segments,  the  last  one  some¬ 
what  overshooting  Thule.  Table  5  gives  the  means  T^  from  Figure  4a,  the  standard 
deviations  o\  from  Figure  4b  and  the  parameters  a .  from  Figure  4c.  Since  the  check 
points  are  evenly'spaced  on  the  chart,  each  weight  is  made  inversely  proportional 
o  the  scale  factor,  that  is, 

s.  =  1  +  sin0.  .  (11) 


The  end  points  XQ  and  Xg  are  given  half  this  weight.  For  each  pair  of  points, 


a  . . 
ij 


(o,i  +  a  j) 


,  =  JUl  -  i> 

lj  10 


radians, 
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Table  5.  The  Straightline  Route  on  the  Polar  Stereographic  Chart,  San 
Francisco  to  Thule,  Winter,  50  mb.  (Segment  divisions  are  0.  1  earth 
radian  in  length.  For  each  checkpoint  the  table  shows  latitude  <f>^,  in  °N; 
mean  temperature  Tj,  in  °C;  standard  deviation  0^,  in  °C;  weighting 
factor  si,  which  is  inversely  proportional  to  the  scale  factor;  and 
parameter**!.  The  terms**  p^,  and  S;S-0.0;,  in  that  order,  are  listed 
in  the  body  of  the  table.  )  J  J  J 


Station  No. 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

latitude 

3H 

42 

47 

52 

57 

62 

66 

71 

75 

78 

Moan  Tj 

SD  a. 

-GO 

-59 

-58 

-57 

-57 

-57 

-57 

-58 

-62 

-64 

3.  3 

3.  7 

4.  2 

5.  0 

6.  2 

7.  4 

8.  5 

10  0 

10.  5 

9.  8 

•I 

0.  H 

1.  7 

1.  7 

1.  8 

1.  8 

1. 9 

1.9 

1.  9 

2.  0 

1.0 

2.  5 

2.  5 

2.  5 

2.  4 

2.  2 

2.  0 

2.  3 

2.  3 

2.  4 

2.  4 

Station  No. 

m 

2.  5 

2.  4 

2.  4 

2.  3 

2.  3 

2.  4 

2.  4 

2.  4 

0 

1.00 

0.  93 

.  77 

.  58 

.  36 

.  20 

.  05 

-.  03 

-.  05 

-.  05 

6.8 

16.  4 

18.  5 

23.  4 

29.  1 

36.  7 

42.  1 

49.  4 

54.  6 

25.  5 

2.  5 

2  4 

2.  3 

2.  3 

2.4 

2.  4 

2.  4 

2.  4 

1 

1. 00 

.  93 

.  74 

.  61 

.  40 

18 

.‘05 

-.  03 

-.  05 

39.  7 

44.  7 

56.  7 

70.  G 

88.  8 

102.  1 

119.  7 

132.  3 

61.  7 

2.  4 

2.  4 

2  3 

2.  4 

2  4 

2.  4 

2  4 

2 

1. 00 

.  94 

.  78 

.  61 

.  36 

.  ia 

.  05 

-.  03 

50.  4 

63.  9 

79.  5 

100.  1 

115.  0 

134.  9 

149.  1 

69.  6 

2.  3 

2  2 

2.  4 

2.  4 

2.  4 

2.  4 

3 

t.  00 

.95 

.  82 

.  58 

.  36 

18 

.  05 

81. 8 

100.  8 

126.  9 

145.  8 

171. 0 

189.  0 

88.  2 

2.  1 

2.  2 

2.  2 

2.  3 

2.  3 

4 

l.  00 

.  95 

.  82 

.  64 

40 

.  21 

125.  4 

157.  9 

181. 4 

212.  8 

235.  2 

109.  8 

2.  2 

2.  2 

2.  2 

2.  3 

5 

1.  00 

.  95 

82 

.  63 

43 

198.  8 

228.  4 

267.  9 

296.  1 

1  38.  1 

2.  3 

2.  4 

2  4 

(i 

1.  00 

.  95 

78 

.  58 

262.  4 

307.  8 

340.  2  1 

158.  8 

2.  4 

2.  4 

7 

1.  00 

.  94 

.  79 

361.  0 

399.  0 

186.  2 

2.  4 

i! 

1.  00 

.  94 

441  0 

205  8 

9 

4 

l .  no 

96.  0 
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from  which  it  is  possible  to  obtain  0“-  from  Figure  3  or  Eq.  (4).  In  Table  5  the 
values  of  are  shown  together  with  the  associated  values  of  SjS.aja. .  Finally, 
the  route  variance  is  obtained  by  cumulative  multiplication: 


r  9 


i=o 


2  2 
'i  ai 


+  2 


E 

i'j 


(12) 


Table  2  shows  the  standard  deviations  estimated  according  to  the  model  (Column  A). 
The  root -mean-square  difference  between  these  estimates  and  estimates  of  the 
FAA  is  0.  4^C,  the  same  as  the  difference  between  the  estimates  of  Columns  B  and  C, 
the  latter  two  being  equally  valid. 

To  test  the  efficiency  of  the  assumption  that  the  route  average  is  normally 
distributed,  the  1 -percentiles  and  5-percentiles  were  computed  by  the  formula 

T*  s  T  +  yo,  (13) 

where  y  *  2.  33  for  the  1 -percentile,  y  8  1,  645  for  the  5-percentile.  This  estimate 
was  compared  with  the  estimates  given  by  the  distributions  in  the  tables  of  the  FAA 
*  studies.  (The  latter  distributions  are  actual  distributions  independent  of  assumptions 
about  their  nature. )  The  results  are  shown  in  Table  6.  The  root -mean-square 
difference  is  1. 1°C  for  the  1 -percentile  and  0. 1®C  for  the  5-percentile. 


Table  6.  Comparison  of  Two  Estimates  of  the  1 -Percentile 
ard  5-Percentile  of  Average  Route  Temperature 


1  -  Percentile 

5-Percentile 

oeason 

nOUIe 

A* 

B+ 

A* 

Bt 

Summer 

SF  —  NY 

SF  —  Thule 

-51.  2 
-45.  2 

-52.  1 
-45.  0 

-51.  9 
-46.  0 

-52.  5 
-45.  6 

Winter 

SF  —  NY 

SF  —  Thule 

-52.  7 
-46.  0 

-53.  9 
-47.  6 

-54.  5 
-49.  8 

-55.  3 
-49.  5 

*by  model 

ffrom  USWB-FAA  studies 


6.  SUMMARY 

To  obtain  the  frequency  distribution  of  the  route  average  temperature  on  a 
constant -pressure  surface,  it  is  assumed  that  the  route  average  has  a  normal 
gaussian  distribution.  The  mean  of  the  route  average  is  given  by: 
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[Eq.  (2)] 


where  there  are  (n+  1)  checkpoints  along  the  route, including  the  point  of  departure 
and  destination.  The  mean  temperature  at  each  checkpoint  is  ar.d  the  weight 
given  to  it,  ,  is  made  directly  proportional  to  the  fraction  of  the  route  over 
which  is  a  representative  temperature.  In  this  paper  (see  page  12), 


j.  =  1  +  sin<^ 


(Eq.  (11)] 


where  0.  is  the  latitude  of  the  ith  checkpoint.  The  variance  of  the  route  average 
is  given  by  [see  Eq.  (12)] 


n,  n  j  n  2 

-2  •  £  si  ^/  £  si) 

1,  j  =  o  /  '  1=0  t 


where  <7.  is  the  standard  deviation  at  the  ith  checkpoint,  and  ^  is  the  correlation 
coefficient  given  by 


fry  *  e*P[-  Vij)  • 


[Eq.  (4)) 


The  distance  y..  in  this  paper  is  the  chart  distance,  in  earth  radians,  given  by 


(Xj  -  X.)2  +  (y.  -  yj)2 


[Eq.  (5)1 


where 


0. 

2cosA{  tan(J  -  . 


0. 

y.  =  2  sin  A.  tan(J  -  ~) 


where  (0.,  A.)  are  the  latitude  and  longitude  of  the  ith  checkpoint.  The  parameter 
a  .j  is  an  average,  that  is, 
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«ij  =  -*2  ^Eq-  (6)1 

where  is  the  parameter  value  at  the  ith  station,  depicted  in  the  (c)  charts  of 
Figures  4  to  11, 


7.  CONCLUSIONS 

The  sample  application  shows  that  the  model,  consisting  primarily  of  a 
parameter  to  obtain  spatial  correlation  of  temperature  on  a  constant -pressure 
surface,  is  an  efficient  tool  for  estimating  the  standard  deviation  of  the  route 
average  temperature  and  an  effective  one  for  estimating  the  5-  or  1 -percentiles. 
There  are  refinements  that  could  have  been  tried  but  would  have  been  unnecessary. 
For  example,  the  correlation  coefficient  might  decay  more  rapidly  in  a  N-S  direction 
than  in  an  E-W  direction,  hence  making  it  desirable  for  the  value  of  the  parameter 
a  to  be  a  vector  combination  of  two  other  values  (say  o  ,  «w>  .  Since  the  simpler 
model  with  only  a  single  parameter  per  station  works  so  well,  attempts  at  refine¬ 
ments  were  abandoned. 


Figure  4a.  Isotherms  of  Mean  Temperature  T.,  50  mb,  December,  January,  February 


Figure  4&.  Isopleths  of  Standard  Deviation  of  the  Temperature  at  50  mb,  Decem¬ 
ber,  January,  February 


I 

i 

* 

} 

i 


Figure  4c.  Isopleths  of  the  Model  Parameter  cr^,  50  mb,  December,  January,  February 
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Figure  5a.  Isotherms  of  Mean  Temperature  T.,  50  mb,  March,  April,  May 


Figure  6a.  Isotherms  of  Mean  Temperature  T.,  50  mb,  June,  July,  August 
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Figure  6c.  Isopleths  of  the  Model  Parameter**.,  50  mb,  June,  July,  August 
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Figure  7q.  Isopleths  of  Model  Parameter  c^,  50mb,  September,  October,  November 
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Figure  8a.  Isotherms  of  Mean  Temperature  T^,  30  mb,  December,  January,  February 

i 


i 
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Figure  8g.  Isopleths  of  Model  Parameter  cr,  30  mb,  December,  January,  February 
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! 


33 
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Figure  10^.  Isopleths  of  Standard  Deviation  cr.,  30  mb,  June,  July,  August 


I 
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Figure  11a.  Isotherms  of  MeanTemperature  T.,  30  mb,  September,  October,  November 


Figure  life.  Isopleths  of  Standard  Deviation  a j,  30 mb,  September,  October,  November 


Figure  11c.  Isopleths  of  Model  Parameter  a.,  30  mb,  September,  October,  November 


Figure  12.  Isopleths  of  the  Difference  D_Between_Mean  Temperatures,  30  mb. 
Winter,  at  Times  0000Z  and  1200Z  (D  =  T12Q0  - 


Figure  13.  Isopleths  of  the  Difference  D  Between  Mean  Temperatures,  30  mb. 
Summer,  at  Times  0000Z  and  1200Z  (D  =  T^qq  "  ^nnoO^ 
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